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Sensitivity  of  the  deconvolution  of  acoustic  transients  to  Green’s 
function  mismatch 


Michael  K.  Broadhead,  Robert  L.  Field,  and  James  H.  Leclere 

Naval  Research  Laboratory,  Stennis  Space  Center.  Mississippi  39529 

(Received  13  July  1992,  revised  22  March  1993;  accepted  22  April  1993) 

In  this  paper  various  measurements  of  acoustic  time  series  recorded  in  the  Atlantic  over 
source-to-receiver  ranges  of  600  to  12  900  m  are  analyzed.  The  transmitted  source  signature  with 
a  monitor  hydrophone  mounted  on  the  source  array  was  also  recorded.  The  signature  distortion 
introduced  by  propagation  effects  was  treated  by  the  use  of  single-channel  deconvolution.  In 
situations  where  the  Green’s  function  structure  is  simple  (e.g.,  direct  arrival  and  surface 
reflection),  single-channel  deconvolution  gave  satisfactory  results.  When  multipath  effects  (due 
to  interaction  with  layered  bottom  sediments)  were  present,  it  was  difficult  to  get  a  good  source 
estimate.  A  way  was  developed  of  perturbing  the  Green’s  functions  such  that  the  source 
estimates  were  guaranteed  to  improve.  In  most  cases  it  was  found  that  very  small  changes  could 
produce  significant  improvement  in  the  source  estimates.  This  sensitivity  was  quantified  by  using 
the  correlation  coefficient.  This  sensitivity  is  not  to  be  confused  with  the  well-known  fact  that  the 
single-channel  deconvolution  problem  is  ill-posed.  That  issue  was  treated  separately. 

PACS  numbers:  43.30. Wi 


INTRODUCTION 

Acoustic  transients  are  signals  that  typically  have 
broadband  spectra  and  a  well-defined  time  du'-ation,  as 
opposed  to,  say,  the  persistent,  narrow-band  signature  of  a 
propeller  blade.  Sources  of  underwater  acoustic  transients 
include  man-made  devices,  sounds  produced  by  natural 
phenomena  such  as  ice  shifting  or  breaking,  and  biologic 
sounds  such  as  those  produced  by  whales  and  other  sea  life. 

The  problem  of  determining  the  nature  of  a  transient 
sound  source,  or  the  classification  problem,  is  currently  of 
interest  to  the  underwater  sound  community.  Much  of  the 
effort  to  date  has  centered  on  the  problem  of  signals  that 
have  not  been  distorted  by  propagation  through  the  envi¬ 
ronment.  Propagation  effects  (or  environmental  distor¬ 
tion)  include  attenuation  and  multipath  effects.  In  this 
study,  we  will  confine  ourselves  to  low  frequency  (<150 
Hz)  phenomena,  so  we  will  expect  most  of  the  attenuation 
to  occur  due  to  bottom  interaction.  Since  we  are  dealing 
with  relatively  short  ranges,  we  will  mainly  be  concerned 
with  multipath  effects.  One  approach  to  classification  in 
the  presence  of  environmental  distortion  effects  is  to  first 
calculate  the  propagation  effects,  and  then  remove  them  in 
a  signal  processing  step.  The  output  would  then  be  an  es¬ 
timate  of  the  transient  source  signature  that  could  be  used 
as  input  into  any  one  of  a  number  of  classification  (pattern 
recognition)  schemes.  We  can  account  for  propagation  ef¬ 
fects  by  computing  the  Green’s  function  G(r,t;r',r').  For  a 
point  source  at  r'  and  a  receiver  at  r.  the  processing  step 
requires  deconvolving  the  pressure  time  series  i/'(r,/)  with 
G{T,ty,t')  for  the  desired  source  signature  S{t).  This  in¬ 
volves  solving  the  following  Fredholm  integral  equation  of 
the  first  kind: 


G(T:r',t~t')S(l’)dt', 


(1) 


where  r  and  r'  are  fixed  for  a  given  case.  It  is  well  known 
that  this  problem  (single-channel  deconvolution)  is 
ill-posed.'"^  This  is  a  different  type  of  sensitivity  from  the 
kind  we  are  studying,  and  we  will  return  to  this  issue  later. 

A  reasonable  question  to  ask  is,  how  accurately  does 
the  Green’s  function  have  to  be  computed  for  the  decon¬ 
volution  to  be  successful?  The  information  that  determines 
the  Green’s  function  includes  source  and  receiver  loca¬ 
tions,  the  sound  speed  and  density  profiles  of  the  water 
layer,  bottom  topography,  and  the  sediment  layer  param¬ 
eters  (velocity,  densit and  attenuation  profiles).  If  these 
environmental  parameters  are  not  known  to  sufficient  ac¬ 
curacy  and  spatial  resolution,  there  will  be  errors  (mis¬ 
match)  in  the  Green’s  function.  The  practical  issue  then 
becomes,  how  sensitively  does  the  source  signature  decon¬ 
volution  problem  (actually,  its  well-posed  extension)  de¬ 
pend  on  G(r,t;r',/')? 

The  work  reported  here  contributes  to  the  study  of 
these  issues  by  using  high  SNR  experimental  data  in  mod¬ 
erate  water  depth  (915  m).  for  the  low-frequency  case,  in 
a  mid  to  short  range  scenario  (600  m- 13  km),  in  which  the 
subbottom  sediment  layers  played  an  important  role  in  the 
multipath  structure.  We  performed  a  sensitivity  analysis  on 
the  deconvolution  of  the  data  with  computed  Green’s  func¬ 
tions  using  the  best  available  geoacoustic  information.  This 
study  was  part  of  a  larger  effort  to  analyze  these  data  while 
simultaneously  evaluating  a  (then)  new  propagation 
model.  ^ 

We  proceeded  with  the  sensitivity  analysis  by  decon¬ 
volving  the  data  with  computed  Green’s  functions  and 
then  perturbing  the  Green’s  functions  in  such  a  way  that 
the  source  estimates  would  be  guaranteed  to  improve.  We 
then  used  the  correlation  coefficient  measure  to  show  that 
significant  increases  in  deconvolution  quality  were  gained 
from  small  changes  (U'-ually)  lo  the  calcidated  Green’s 
function.  This  leads  us  to  su.spect  that  the  gcoacoustic  pa- 
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FIG.  1.  Experiment  geometry.  Location — Blake  Plateau;  water  a.pth — 
915  m;  source — HLF-2AH  towed  source;  bandwidth — 25*  to  150-Hz  lin¬ 
ear  sweep;  array — 15-element  vertical  array. 
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rameters  and  source  location  may  need  to  be  known  to  a 
greater  accuracy  than  will  be  feasible,  at  least  in  environ¬ 
ments  where  the  multip’th  structure  becomes  sufficiently 
complex. 

In  the  rest  of  this  paper  we  will  describe  the  data  set 
we  have  worked  with,  the  modeling  that  was  done,  and  the 
initial  deconvolution  results  using  the  model  Green’s  func¬ 
tion.  We  will  then  present  the  optimization  method  that 
will  allow  us  to  perturb  the  model  Green’s  function  to  get 
better  deconvolution  results,  and  show  results  from  calcu¬ 
lations  with  the  data.  Finally  we  will  summarize  our  results 
and  present  some  conclusions. 

I.  EXPERIMENT  DESCRIPTION 

The  experiment  from  which  our  data  were  obtained 
was  conducted  in  the  vicinity  of  27.5  °N,  78.3  ’W  (Blake 
Plateau)  in  the  Atlantic  Ocean  (see  Fig.  1).  The  water 
bottom  was  fairly  smooth  and  the  water  depth  was  approx¬ 
imately  915  m.  The  experiment  included  measurements  us¬ 
ing  a  23-  to  150-Hz  linear  sweep  pulse  produced  by  an 
HLF-2AH  towed  acoustic  source.  The  signals  were  re¬ 
corded  on  a  stationary  15-element  vertical  line  array  with 
9-m  hydrophone  spacing.  The  top  hydrophone  was  at  a 
depth  of  Z=129  m  and  the  bottom  hydrophone  was  at 
Z=250  m.  Refer  to  Field  and  Leclere^  for  more  detail. 

II.  DATA  PREPARATION  AND  MODELING 

In  Fig.  2  we  display  the  (correlated)  time  series  for  the 
five  source-to-receiver  ranges  (R),  600,  1500,  43(X),  7900, 
and  12  900  m.  The  hydrophone  depth  is  Z=  129  m.  These 
data  have  been  cross  correlated  with  the  measured  source 
signature  from  a  monitor  hydrophone  to  improve  the  SNR. 


FIG.  2.  Measured  transmissions  at  receiver  depth  Z=129  m  for  the  five 
ranges  considered.  Data  have  been  correlated  with  the  measured  source 
signature. 

In  the  R  =  600  m  and  /f=  1500  m  cases,  the  arrival  struc¬ 
ture  is  primarily  composed  of  direct  arrival  and  surface 
bounce.  There  are  some  low  amplitude  bottom  interactions 
past  0.8  s.  However,  in  the  other  three  cases,  /?  =  4.3,  7.9, 
and  12.9  km,  there  is  significant  bottom  interaction 
present,  which  plays  an  important  role  in  this  study  ( refer 
to  Field  and  Leclere^).  In  Fig.  3  we  show  the  (correlated) 
measured  source  signature.  This  quantity  plays  the  role  of 
S{t)  in  our  calculations.  It  is  also  the  ground  truth  or 
standard  by  which  our  source  estimates  will  be  judged.  In 
Fig.  4,  we  display  the  computed  model  Green’s  functions 


FtG.  3.  The  source  signature  (after  correlation)  that  was  transmitted 
(measured  from  a  monitor  hydrophone). 
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FIG.  4.  Computed  (TDPE)  band-limited  Green’s  functions  for  the 
129-m  hydrophone. 

for  the  corresponding  measurements.  The  model  used  was 
the  time  domain  parabolic  equation  (TDPE)  model,'*’’ 
Bathythermographs  were  used  to  obtain  water  column 
sound  speeds  while  sediment  velocities  and  densities  were 
extrapolated  from  data  gathered  at  Deep  Sea  Drilling 
Project  sites  near  the  experiment.  These  environmental  pa¬ 
rameters  were  input  into  TDPE  along  with  a  desired 
source-to-receiver  range,  producing  a  band-limited  (0  to 
150-Hz)  Green’s  function  as  output.  Refer  to  Field  and 
Leclere’  for  more  detail. 

In  Fig.  5  we  show  the  simulations  produced  by  con¬ 
volving  the  source  S(t)  in  Fig.  3  with  the  Green’s  func¬ 
tions  in  Fig.  4.  The  simulations  are  overlaid  (dashed)  with 
the  measurements  (solid)  in  Fig.  2.  For  a  q’uantitative 
comparison,  note  the  correlation  coefficients  (y)  that  are 
provided,  y  is  defined  in  the  usual  way  .so  that 
where 


y=  X  (2) 

i  '  *  /  / 

and  are  means.  Notice  that  for  /?  =  600  m,  0.883 
and  for  i?  =  1500  m,  y=0.827.  This  demonstrates  that  the 
model,  environmental  input,  and  the  various  assumptions 
made  are  ail  reasonable  for  at  least  some  cases.  However, 
for  /?>4.3  km,  the  simulation  does  not  match  the  data  as 
well  (y<0.638).  This  is  presumed  to  be  at  least  partly  due 
to  the  fact  that  significant  bottom  interaction  is  present. 
The  bottom  interaction  presents  a  problem  because  we 
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FIG.  5.  Measurements  in  Fig.  2  (solid)  compared  with  simulatirtis 
(dashed)  produced  by  convolving  source  signature  in  Fig.  3  with  Green’s 
functions  in  Fig.  4.  The  correlation  coefficient  provides  a  quantitative 
comparison. 

have  limited  information  on  the  sediment  layer  geophysical 
parameters. 


III.  DECONVOLUTION  ALGORITHM  AND  RESULTS 

As  mentioned  earlier,  the  importance  of  the  Green's 
function  from  the  point  of  view  of  the  transient  classifica¬ 
tion  problem  is  that  it  gives  us,  along  with  the  received 
signal,  enough  information  to  solve  Eq.  ( 1 )  for  the  source 
signature  S(t).  We  will  present  some  deconvolution  re¬ 
sults,  but  first  we  will  discuss  some  of  the  theoretical  and 
computational  aspects.  We  have  already  pointed  out  that 
the  deconvolution  problem  represented  by  Eq.  ( 1 )  is  math¬ 
ematically  ill-posed.  Exactly  what  does  this  mean?  Accord¬ 
ing  to  Tikhonov  and  Arsenin’  and  Tarantola:^ 

Definition  1:  If  a  solution  S[t)  exists,  is  unique,  and  is 
a  continuous  function  of  the  initial  data  then  the 

problem  is  well-posed.  Otherwise,  it  is  ill-posed. 

Definition  2:  Small  changes  in  the  initial  data 
yield  arbitrarily  large  changes  in  the  solution  S(/). 

Tikhonov  and  Arsenin’  point  out  that,  in  practice, 
many  mathematical  problems  of  physical  interest  are  ill- 
posed.  They  and  others  have  published  techniques  for  pro¬ 
viding  well-posed  extensions  of  these  problems.  Tikhonov 
regularization,'  ^  pre-whitening,  and  singular  value  de¬ 
composition  (SVD)’’’"’  are  some  examples.  Before  we  pro¬ 
ceed  further  with  this  discussion,  it  will  be  useful  for  us  to 
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recast  Eq.  ( 1 )  in  a  discrete  form  that  is  amenable  to  com¬ 
puter  implementation. 

Let  R  and  Z  be  fixed.  Then  the  discrete  (computed) 
Green’s  function  can  be  written  as  a  column  vector 
where  the  subscript  indexes  time  sam¬ 
ples.  We  can  also  write  this  discrete  Green’s  function  in  a 
convolutional  matrix  form," 
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If  Sj  is  the  column  vector  for  the  unknown  source 
estimate  and  d  is  the  data  vector  for  our  choice  of  R  and  Z, 
then  the  discrete  form  of  the  convolution  equation  can  be 
written, 

Gs.=d.  (4) 

As  is  well  known,  this  linear  system  is  over  determined  and 
may  have  no  solution."  A  common  approach  to  construct¬ 
ing  a  solution  is  to  use  the  method  of  least  squares.  We 
define  an  error  vector  c  as 

e=Gs^-d.  (5) 

We  then  seek  so  that  the  square  length  is  a  mini¬ 
mum.  The  normal  equativ  ns  obtained  from  setting  the  gra¬ 
dient  to  zero  are 

G^Gs*=G^d.  (6) 

This  linear  system  is  evenly  determined  but  may  be  ill- 
conditioned.  This  is  the  matrix  version  of  the  fact  that  the 
original  integral  equation  was  ill-posed.  In  order  to  provide 
a  well-posed  extension,  we  chose  the  method  of  pre- 
whitening,  which  has  been  shown’  to  be  a  form  of 
Tikhonov  regularization.'*  We  initially  used  SVD,  but 
found  the  extra  computational  burden  unnecessary.  The 
form  of  the  normal  equations  with  pre-whitening  are 

(G’'G+AI)s=G’’d,  (7) 

where  A=e<^o,  0o  Ihe  main  diagonal  element  (they  are 
all  the  same),  e  is  a  small  positive  number,  and  I  is  the 
identity  matrix.  The  idea  behind  pre-whitening  is  that  it 
prevents  any  eigenvalues  of  the  coefficient  matrix  from  be¬ 
coming  too  small  by  giving  them  all  a  positive  shift.  Here, 
A  can  also  be  interpreted  as  a  Lagrange  multiplier  for  the 
addition  of  a  finite  filter  energy  constraint.’*  Since  the  nor¬ 
mal  equations  have  a  Toeplitz  structure,"  we  were  able  to 
use  the  efficient  Levinson  recursion  method  for 
solution."’"  As  a  further  processing  step,  we  applied  a 
hi-cut  filter  to  preserve  only  the  original  Green’s  function 
bandwidth. 

In  Fig.  6  we  have  displayed  the  results  of  deconvolving 
the  five  measurements  in  Fig.  2  with  the  Green’s  functions 
in  Fig.  4.  The  results  are  the  dashed  plots  and  are  overlaid 
with  the  measured  source  signature  (solid)  that  was  pre¬ 
viously  shown  in  Fig.  3.  We  have  provided  the  correlation 
coefficients  for  a  quantitative  comparison.  Note  that  they 
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FIO.  6.  Source  signature  (solid)  compared  with  source  estimate 
(dashed)  produced  by  deconvolving  measurements  in  Fig.  2  with  Green's 
functions  in  Fig.  4. 

are  comparable  to  the  correlation  coefficients  in  Fig.  5, 
where  we  compare  simulations  and  measurements.  We  will 
see  the  close  connection  between  these  two  quantities  again 
in  this  paper.  It  is  not  suprising  that  the  match  between  the 
simulations  S{t)*G(t)  (*  means  convolution)  and  the 
data  are  a  good  predictor  of  the  quality  of  the  deconvolu¬ 
tions.  We  should  point  out  that  the  source  estimates  for 
/?  =  600  m  and  /?  =  1500  m  are  probably  of  a  “usable” 
quality  in  a  classification  sense;  /I  =  12.9  km  may  be  mar¬ 
ginally  usable;  /?=4300  and  7900  m  are  probably  not  us¬ 
able.  So,  a  simple  fall  off  in  quality  with  range  is  not  a 
justified  conclusion. 

IV.  MISMATCH  PROBLEM  AND  SOLUTION 

It  was  of  interest  for  us  to  determine  how  much  the 
environmental  parameters  would  have  to  change  to  pro¬ 
duce  correlation  coefficients  that  were  comparable  to  the 
R  ~  600  m  and  R  =  1 500  m  cases.  To  solve  this  problem 
directly  was  computationally  prohibitive  in  terms  of  the 
resources  available.  As  an  alternative,  we  chose  to  perturb 
only  the  Green’s  function  itself,  instead  of  the  environmen¬ 
tal  parameters  the  Green’s  function  is  calculated  from. 

In  order  to  perturb  the  Green’s  function  in  a  meaning¬ 
ful  way,  we  set  up  a  double  objective  optimization  problem, 
i.e.,  one  with  two  goals.  If  our  current  G(t)  (where  we 
have  fixed  r  and  r')  does  not  produce  a  simulation 
S{t)*G(t)  with  good  correlation  to  the  measured  data, 
then  we  seek  a  new  G'(t)  that  does.  However,  we  do  not 
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want  G'(t)  to  be  too  different  from  G{t),  which  is  our  tie 
to  the  physics  of  the  problem.  In  other  words,  we  do  not 
want  G'(r)  to  be  arbitrary.  The  problem  was  set  up  in  the 
following  way;  for  some  fixed  R  and  Z,  let  the  column 
vector  g  be  the  perturbed  Green’s  function  we  are  seeking 
and  let  be  the  model  calculated  Green’s  function.  Let  S 
be  the  matrix  convolutional  form  representing  the  discrete 
source  signature.  Recall  that  d  is  the  correlated  measure¬ 
ments  vector  at  range  R  and  hydrophone  depth  Z.  Then 
the  mismatch  (oi  difference)  between  the  perturbed  simu¬ 
lation  and  the  measured  data  will  be 

e,=Sg-d.  (8) 

The  mismatch  between  the  perturbed  Green’s  function  and 
the  model  Green’s  function  will  be 


^2=S-Sm  (9) 

Since  we  wish  to  see  how  much  it  is  possible  to  reduce 
the  mismatch,  we  need  a  quantity  that  measures  the  sizes 
of  these  vectors.  Again  we  choose  the  least  squares  criteria. 
We  seek  a  g  to  minimize  the  functional: 

E{g)=ejex+0^e^ei.  (10) 

According  to  Scales,''*  a  quadratic  form  can  be  written  as, 

/’(x)=5X^Ax+b^x-|-c,  (11) 

and  its  gradient  can  be  written, 

VF(x)=Ax-(-b.  (12) 

If  we  substitute  Eqs.  (8)  and  (9)  into  Eq.  (10),  and  rear¬ 
range  terms  we  get 

£(g)  =kl^[2S^S-)-2/82i]g+ [  _s^d-^'g^]  ^g+c. 

(13) 

Comparing  (13)  to  (11),  and  then  using  ( 12)  we  can  write 

VF(g)=2[S^S+i3'l]g-f  [-S^d-/?'g^].  (14) 

Upon  setting  the  gradient  to  zero,  we  have  the  normal 
equations; 

(S’^S+R^Dg^HS^d+Phu)-  (15) 

This  system  was  also  solved  using  Levinson  recursion. 
Comparison  of  Eq.  (15)  with  Eq.  (7)  shows  how  0^  plays 
a  role  similar  to  A.  Although  we  also  used  pre-whitening  in 
our  numerical  solution  of  (15),  it  turned  out  to  be  needed 
only  when  0^  was  zero  (or  sufficiently  small). 

Let  us  briefly  examine  what  the  functional  in  Eq.  ( 10) 
means.  If  we  consider  just  the  first  term  alone,  then  the 
problem  becomes  one  of  finding  a  g  that,  undrr  convolu¬ 
tion  with  the  source,  produces  a  good  match  to  the  mea¬ 
surements  (in  a  least-squares  sense).  If  we  consider  Just  the 
second  term  alone,  the  problem  is  that  of  finding  a  g  that  is 
most  similar  to  g^.  This  problem  has  the  trivial  solution 
g=gAf-  When  we  put  the  two  terms  together,  however,  we 
find  we  have  a  double  objective  optimization  problem.  If 
the  goals  are  conflicting,  the  solution  will  be  a  compromise. 
The  term  that  receives  the  most  emphasis  in  the  compro¬ 
mise  is  governed  by  the  trade-off  parameter  0^.  As  0  ap¬ 
proaches  infinity,  g  approaches  g^j, .  As  ^  approaches  0,  g 


R  =  600  m,  Z  =  129  m 


R=  1500  m,Z=  129  m 


R  =  4300  m,  2  =  129  m 


R  =  7900  m,  Z  =  129  m 


FIG.  7.  Beta  versus  correlation  coefficients  between  the  TDPE  mode! 
Green’s  functions  and  the  optimized  Green's  functions  (solid);  the  model 
simulations  and  the  optimized  simulations  (dashed);  and  the  source  sig¬ 
nature  and  the  optimized  source  estimates  from  deconvolution  (long 
dash) 


approaches  the  unconstrained  deconvolution  solution.  An 
other  way  to  view  Eq.  ( 10)  is  that  we  are  adding  a  con 
straint  that  g  be  close  to  . 


V.  OPTIMIZATION  RESULTS 

In  Fig.  7,  we  show  a  summary  of  results  from  the 
double  objective  problem.  For  a  particular  R  and  Z,  and  a 
particular  0  value,  the  quantity  g  was  produced  by  solving 
the  linear  system  in  Eq.  (15).  This  “optimized”  Green’s 
function  was  then  convolved  with  s  to  produce  a  simula¬ 
tion  of  the  correlated  data  and  was  used  to  deconvolve  d  to 
produce  a  source  estimate  Sp.  The  correlation  coefficients 
between  the  model  and  optimized  Green’s  function 
(solid),  the  measured  data  and  the  optimized  simulation  of 
the  data  (dashed),  and  the  source  signature  and  the  opti¬ 
mized  source  estimate  from  deconvolution  (long  dash)  are 
then  plotted  against  0.  These  curves,  for  each  of  the  five 
ranges  considered,  are  shown  in  Fig.  7.  We  can  immedi¬ 
ately  see  several  things.  First,  as  noted  earlier,  simulation 
quality  and  deconvolution  quality  generally  follow  one  an¬ 
other  closely.  The  exception  is  R~7900,  for  40</8<120. 
Second,  as  expected,  the  simulation  and  deconvolution 
quality  increase  as  g  is  allowed  to  depart  from  g^.  This 
effect  is  more  pronounced  for  the  last  three  ranges  than  for 
the  first  two.  For  0>3OO,  0  is  effectively  cc,  and  g  =  g,v 
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FIG.  8.  For/9=80,  /f=7900  m,  and  Z=  129  m  we  have:  (a)  comparison 
of  model  Green’s  function  (solid)  and  optimized  Green’s  function 
(dashed),  /=0.987;  (b)  comparison  of  measured  data  with  optimized 
simulation,  ^=0,551;  and  (c)  comparison  of  source  signature  with  source 
estimate  obtained  by  deconvolution  of  measurement  with  the  optimized 
Green’s  function,  >'=0,716. 

We  take  particular  note  of  the  crossover  point  in  the 
curves,  which  occurs  at  an  average  value  of  P^29.  We 
have  essentially  already  seen  the  results  for  /3>300,  so  we 
will  pick  three  interesting  P  values  for  a  particular  range 
and  look  at  them  in  more  detail.  Because  the  most  im¬ 
provement  was  gained  for  R  =  790Q  m,  we  will  use  this 
range.  For  reasons  we  will  explain,  we  will  display  the 
results  for  i3=80,  P=29,  P~0. 

Let  us  now  consider  Fig.  8.  In  Fig,  8(a),  we  display 
the  model  (solid)  and  optimized  (dashed)  Green’s  func¬ 
tions  for  i?=7900  m,  Z=  129  m,  and  a  value  of  0=SO. 
Notice  that  the  departure  of  g  from  g^,  is  very  slight.  This 
is  quantified  by  the  correlation  coefficient  value  of 
7=0.987.  The  simulation  produced  by  convolving  g  with  s 
is  shown  in  Fig.  8(b).  Here,  the  correlation  has  improved 
from  a  value  of  7=0.314  (refer  to  Fig.  5)  to  a  value  of 
7=0.551.  Finally,  the  deconvolution  produced  source  es¬ 
timate  is  shown  in  Fig.  8(c)  [overlaid  with  the  source  sig¬ 
nature  (solid)].  Compare  this  to  the  /?  =  7900  m  case  in 


(a) 


TIME  is) 

FIG.  9.  For/7=29,  7?  =  7900  m,  and  2=129  m  we  have:  (a)  comparison 
of  model  Green’s  function  (solid)  and  optimized  Green’s  function 
(dashed),  >'=0.908;  (b)  comparison  of  measured  data  with  optimized 
simulation,  y=0.877;  and  (c)  comparison  of  source  signature  with  source 
estimate  obtained  by  deconvolution  of  measurement  with  the  optimized 
Green’s  function,  }'=0.935. 

Fig.  6.  The  source  estimate  quality  has  increased  from  a 
value  of  7=  0.364  to  a  value  of  7=0.716.  The  reason  we 
displayed  the  /?=80  case  was  that  we  have  somewhat  sub¬ 
jectively  picked  7=0.7  as  a  cutoff  point  in  a  source  esti¬ 
mate  being  usable  in  a  classification  sense.  Thus,  0=80 
allows  us  to  investigate  the  marginal  case. 

The  significance  of  0=29  is,  we  recall,  that  it  is  the 
crossover  value.  This  means  that  about  equal  weight  is 
being  given  to  both  terms  in  the  optimization  problem.  In 
Fig.  9(a),  we  show  the  g  vs  gv>,  with  7  =  0.908.  Here,  the 
departure  is  more  pronounced,  but  the  two  quantities  are 
still  very  similar.  Notice  that  some  of  the  amplitudes  and 
delays  of  the  impulses  have  changed  slightly.  This  is  the 
sort  of  perturbation  we  would  expect  for  a  slight  change  in 
geometry  {R  and  Z)  or  environmental  parameters  (e.g., 
water  column  or  bottom  sediment  layer  sound  speeds).  In 
Fig.  9(b)  we  show  the  simulation  comparison.  Here, 
7=0.877,  which  should  lead  us  to  expect  a  good  deconvo¬ 
lution.  Indeed,  in  Fig.  9(c),  we  see  that  the  source  estimate 
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FJG.  10.  Fori?=0,  /?  =  7900  m,  and  Z=129  m  we  have;  (a)  comparison 
of  model  Green’s  function  (solid)  and  optimized  Green’s  function 
(dashed)  y=0.008:  (b)  comparison  of  measured  data  with  optimized 
simulation,  y=0.997;  and  (c)  comparison  of  source  signature  with  source 
estimate  obtained  by  deconvolution  of  measurement  with  the  optimized 
Green’s  function,  ^=0.996. 


is  quite  good,  with  a  value  of  j'=0.935.  The  significance  of 
P  is  that  it  is  not  the  best  we  can  do  on  the  simulation  and 
deconvolution,  but  that  it  represents  about  the  best  we  can 
do  without  letting  g  change  too  much  with  respect  to  . 
To  demonstrate  this  in  a  more  quantitative  way,  lets  look 
at  another  extreme,  namely,  P=0.  In  Fig.  10(a),  we  show 
the  g  (dashed)  produced  by  the  requirement  that  under 
convolution,  it  should  be  the  best  possible  match  to  the 
data  (i.e.,  P~0).  There  is  virtually  no  correlation  between 
it  and  guf,  i.e.,  7=0.008.  In  Fig.  10(b),  we  see  that  the 
simulation  is  almost  perfect,  with  7=0.997;  and  in  Fig. 
10(c)  we  have  an  almost  perfect  source  estimate,  with 
7  =  0.996. 

Having  established  the  usefulness  of  (3=29,  let  us  now 
go  back  and  look  at  that  case  for  all  five  ranges.  In  Fig.  11, 
we  display  the  model  Green’s  functions  (solid)  versus  the 
optimized  functions  (dashed)  for  the  five  ranges  consid¬ 
ered,  along  with  the  corresponding  correlation  coefficients. 
The  .significant  feature  to  notice  is  that  the  y  values  are  all 
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FIG.  11.  Comparison  of  model  Green's  functions  (solid)  with  optimized 
Green’s  functions  (daslied)  (0  =  29)  for  the  five  ranges  considered. 


quite  high  (in  the  0.9’s).  This  means  we  have  done  very 
little  perturbing.  Now,  in  Fig.  12,  we  display  the  corre¬ 
sponding  simulations  (cf.  Fig.  11).  Notice  that  here  also 
the  7  values  are  quite  good.  Finally,  in  Fig.  1 3,  we  have  the 
source  estimate  (cf  Fig.  6).  The  source  estimates  are  all 
very  good,  with  y  values  mostly  in  the  0.9’s. 

In  Fig.  14  we  summarize  the  foregoing  results  by- 
graphing  the  change  in  correlation  coefficient  between 
P=<x>  and  J3=29  for  the  g’s,  the  simulations,  and  the 
source  estimates  for  each  of  the  five  ranges  considered 
(here,  numbered  1  through  5).  This  figure  summarizes  in  a 
succinct  and  quantitative  manner  the  results  we  have  de¬ 
scribed  above.  For  example,  notice  that  deconvolution 
quality  and  simulation  quality  track  well.  Also,  the  first 
two  ranges  did  not  show  much  change,  but  the  last  three 
did  (being  of  poor  quality  to  begin  with).  The  main  point 
of  this  figure  is  that  the  change  in  y  for  g  with  respect  to 
gjif,  was  never  very  much  (<0.1),  whereas  the  corre¬ 
sponding  Ay's  for  the  source  estimates  and  simulations 
were  considerable  (  ~0.6).  This  has  led  us  to  suggest  that 
there  is  a  .sensitivity  to  the  exact  structure  of  an  impulsive 
type  Green’s  functions  when  significant  multipath  is  in¬ 
volved.  This  idea  would  have  potentially  significant  impact 
for  the  transients  classification  problem. 

Vi.  SUMMARY 

To  summarize  the  evidence  for  sensitivity  in  a  more 
succinct  way,  we  can  display  the  means  and  the  maxima 
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FIG.  12.  Comparison  of  measurements  (solid)  with  optimized  simula¬ 
tions  (dashed)  (/3=29)  for  the  five  ranges  considered. 
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FIG.  13.  Comparison  of  source  signature  (solid)  with  source  estimates 
(dashed)  obtained  from  deconvolving  measurements  with  optimized 
Green’s  functions  (p—29)  for  the  five  ranges  considered. 


RANGE  (no.) 


FIG.  14.  The  change  in  correlation  coefficients  between  ^  =  300  (model) 
and ^=29  for  ( I)  the  Green's  functions  (  4  ),  and  the  simulations  (jc)  for 
the  five  ranges  considered  (numbered  1-5  as  range  increases)  The  means 
are  over  the  last  three  ranges. 

for  the  curves  in  Fig.  14.  The  averages  over  the  last  three 
ranges  are  (Ay)qp=^0M,  <At'>s,m=0.4I,  and 

<Ay)DCN=0-35,  and  the  maxima  can  be  written  Max(Ay) 
=  {0.09  (GF),  0.56  (SIM),  0.57  (DCN)}.  This  indicates 
that  relatively  large  changes  (Ay~0.6)  in  the  simulation 
and  source  estimate  quality  can  be  caused  by  relatively 
small  perturbations  (Ay~0.1)  to  the  corresponding 
Green’s  functions  when  multipath  is  significant. 

VII.  CONCLUSIONS 

We  have  shown  a  sensitivity  of  the  acoustic  transients 
deconvolution  problem  to  small  errors  in  the  Green’s  func¬ 
tion  when  the  multipath  structure  is  sufficiently  complex. 
Our  multipath  is  mainly  due  to  interaction  with  the  sedi¬ 
ment  layers  in  the  bottom,  which  is  where  we  have  the  least 
information.  We  conclude  that,  in  this  kind  of  scenario, 
estimates  of  accuracy  requirements  in  the  environmental 
parameters  (especially  the  subbottom)  and  in  the  source 
location  should  be  an  important  part  of  any  deconvolution 
scheme  involving  multipath  (especially  bottom  interact¬ 
ing).  We  reiterate  that  the  well-known  sensitivity  due  to 
the  ill-posed  nature  of  the  problem  is  a  different  effect.  It 
consists  of  small  changes  to  the  data  d  causing  large 
changes  to  the  source  estimate  s^ .  We  have  dealt  with  this 
problem  by  replacing  ttie  ongmal  ill-posed  problem  with  a 
well-posed  extension  (in  the  Tikhonov  sense).  Further¬ 
more,  the  simulations  (see  Fig.  14)  show  a  similar  sensi¬ 
tivity  to  the  deconvolutions.  The  simulations  were  pro¬ 
duced  by  convolution,  a  well-posed  mathematical 
operation. 
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